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Abstract. We summarize a study where we test the hypothesis that local black holes 
(BH) are relics of AGN activity. We compare the mass function of BHs in the local 
universe with that expected from AGN relics, which are BHs grown entirely with mass 
accretion during AGN phases. The local BH mass function (BHMF) is estimated by 
applying the well-known correlations between BH mass, bulge luminosity and stellar 
velocity dispersion to galaxy luminosity and velocity functions. The density of BHs in 
the local universe is p B H = 4.6±i;| h%, 7 x 10 5 M Mpc~ 3 . The relic BHMF is derived 
from the continuity equation with the only assumption that AGN activity is due to 
accretion onto massive BHs and that merging is not important. We find that the relic 
BHMF at z = is generated mainly at z < 3. Moreover, the BH growth is anti- 
hierarchical in the sense that smaller BHs (Mbh < 10 7 AIq) grow at lower redshifts 
[z < 1) with respect to more massive ones {z ~ 1 — 3). Unlike previous work, we find 
that the BHMF of AGN relics is perfectly consistent with the local BHMF indicating 
the local BHs were mainly grown during AGN activity. This agreement is obtained 
while satisfying, at the same time, the constraints imposed by the X-ray background. 
The comparison with the local BHMF also suggests that the merging process is not 
important in shaping the relic BHMF, at least at low redshifts [z < 3). Our analysis 
thus suggests the following scenario: local BHs grew during AGN phases in which 
accreting matter was converted into radiation with efficiencies e — 0.04 — 0.16 and 
emitted at a fraction A = 0.1 — 1.7 of the Eddington luminosity. The average total 
lifetime of these active phases ranges from ~ 4.5 x 10 8 yr for Mbh < 1O 7 M0 to 
~ 1.5 x 10 s yr for Mbh > 10 9 M Q . 



1 Introduction 



The standard paradigm for Active Galactic Nuclei is that they are powered 
by mass accretion onto a massive BH (Mbh ~ 10 6 — 10 10 M©). Combined 
with the observed evolution of AGN, this implies that many (if not all) nearby 
galaxies should host a BH in their nuclei as relic of past AGN activity. BHs are 
detected in ~ 40 galaxies and their mass correlates with host galaxy structural 
parameters like bulge luminosity/mass [8,11] and stellar velocity dispersion [4,5]. 
An important open question is if local BHs are relics of AGN activity (i.e. grown 
entirely with mass accretion during AGN phases) or if other processes, such as 
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merging, play an important role. This can be answered by comparing the BHMF 
of local BHs with that expected from AGN relics [10,22]. In recent work [22,3] a 
discrepancy in the BHMF at high masses (Mbh > 10 8 M©) has been found: more 
AGN relics are expected than predicted by the local BHMF. This discrepancy 
can be reconciled by assuming accretion efficiencies larger than the canonically 
adopted value of e = 0.1, i.e. e > 0.2. A refinement of the analysis by [22] is also 
presented in [23] who find e > 0.1. High efficiencies are also required from the 
comparison of /?bh derived from the X-ray Background (XRB) and from local 
BHs [2]. Such high efficiencies, if confirmed, would imply that most, if not all 
BHs, should be rapidly rotating. 

In this paper we investigate the possibility that massive black holes in nearby 
galaxies are relics of AGN activity by comparing the local BHMF with that 
of AGN relics. The only assumption is that AGN activity is caused by mass 
accretion onto the central BH. The work is described in detail in [9] and here we 
focus on a few among the more critical and important issues. Similar analysis, 
reaching conclusions analogous to ours are presented in [12,17] and by the same 
authors in these proceedings. 

2 The Mass Function of Local Black Holes 

The mass function of local BHs can be estimated by simply convolving the ex- 
isting galaxy luminosity [</>(£)] or velocity functions [0(c)] with the MeH-ibui 
and Mbh-o* relations respectively. One should apply corrections to convert from 
total to bulge luminosity in the first case and should take into account the in- 
trinsic dispersion (if any) of the MsH-host galaxy relations. In Fig. la we verify 
that the Mbh-c* and Msn-Lbui relations applied to the galaxy luminosity or 
velocity function [1,18] provide the same BHMF within the uncertainties (esti- 
mated with many Montecarlo realizations of the BHMF). Clearly, the necessary 
condition is that the two relations have the same intrinsic dispersion since very 
different BHMFs are derived when MBH-^bui has dispersion 0.5 in log Mbh at 
given Lbui and Mbh-c* has intrinsic dispersion. This confirms the result by 
[11] that all correlations Mbh - host-galaxy-properties are equally good, i.e. they 
have similar intrinsic dispersion. In Fig. lb we plot the estimate of the local BH 
mass function obtained considering galaxies from all morphological types. The 
density in local BHs is p Bli = 4.6(-1.4; +l.9)(h/0.7) 2 x 10 5 MqMpc -3 . We 
have used the galaxy luminosity functions by [7,16,13] and the galaxy velocity 
function by [18]. 

Our estimate of the local BH density is a factor ~ 1.8 larger than the estimate 
by [22] . pbh is increased by ~ 30% when taking into account an intrinsic disper- 
sion for MsH-^bu] and Mbh-c* (0.3 in log Mbh at constant <r* or Lbui)- Pbh is 
also increased by ~ 50% because we have used the zero points of the MBH-ibui 
and Mbh-c* correlations determined by [11]. These are a factor 1.5 larger than 
those used by [22] (see [20]) because they were determined by considering only 
secure BH detections, where the BH sphere of influence is resolved by the obser- 
vations. Fig. 8 of [6] show that measurements obtained with data either resolving 
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Fig. 1. (a) Local BHMF for early type galaxies based on the SDSS sample of [1]. The 
shaded area and error bars ("barred") indicate la uncertainties. The A indicate the 
assumed intrinsic dispersions of the Mbh-c* or AfeH-ibui relations, (b) Best estimate 
of the local BHMF (shaded area) compared with the BHMF of AGN relics obtained 
using the luminosity function by [21], corrected for the missing Compton-thick AGNs. 

or not resolving the BH sphere of influence provide similar Mbh values, albeit 
with much larger errorbars in the latter case. However, from the same figure it 
can be evinced that on average BH measurements where the sphere of influence 
is not resolved are underestimated by a factor ~ 2, and this fully accounts for 
the larger zero points found by [11]. Apart from the larger zero points, by ex- 
cluding the non-secure BH detections [11] find that Mbh-c* and AfeH-^bui have 
the same dispersion which, as we have just shown, is independently confirmed 
by the requirement of obtaining the same BHMF from A/bh-o* and A/BH-^bui 
applied to [</>(cr)] and [</>(£)], respectively. 

3 The Mass Function of AGN Relics 

The Mass Function of AGN relics is estimated with the continuity equation 
which relates the relic BHMF, N(M,t), to the AGN luminosity function (LF), 
4>(L,t) [19]. The only assumption is that AGNs are powered by mass accretion 
onto a massive BH and that we can neglect merging of BHs. The efficiency of 
mass-to-energy conversion is e and the BH is emitting at the fraction A of the 
Eddington luminosity. L in the AGN LF must be the 'bolometric' luminosity. 
To obtain L, we derive and use bolometric corrections which do not take into 
account the IR radiation (reprocessed UV radiation). Thus, they are a factor 
~ 30% lower than the values used by previous authors. We consider the Hard X- 
ray luminosity function by [21], corrected for the missing Compton-thick AGNs 
(factor ~ 1.6) and we apply a bolometric correction to obtain (f>(L,t). Assuming 
that at z = 3 all BHs are active (this initial condition does not affect the final 
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Fig. 2. Locus where e and A provide the 
best match between local and relic BHMFs. 
The solid and dashed lines indicate an av- 
erage deviation of 1 and 0.7a" between the 
BHMFs. The diamond marks the e,A val- 
ues providing the best agreement. 
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results), we can estimate the relic BHMF (with s = 0.1 and A = 1) and compare 
it with the local BHMF (Fig. lb). The local BHMF and the relic BHMF arc in 
good agreement within the uncertainties. Thus, it is unlikely that merging can 
play a major role in shaping the BHMF for z < 3. The [21] LF, corrected for 
the missing Compton-thick AGNs can also reproduce the XRB spectrum and 
source counts, thus satisfying the constraints imposed by the XRB. In particular, 
the disagreement found by [2] between the density of local massive BHs and 
that inferred from the X-ray background light can be reconciled noting that the 
average redshift of the sources making the XRB is not (z) ~ 2 but (z) ~ 1, 
as shown by the redshift evolution of the [21] LF. It is worth stressing the 
importance of the XRB constraint which effectively removes one free parameter 
in this analysis, i.e. the fraction of obscured AGNs. Clearly, a more refined fit of 
the XRB spectrum is required and is the subject of future work. 



4 Accretion efficiency and L/L Edd 

Accretion efficiency e and Eddington ratio A are the only free parameters for 
the relic BHMF and Fig. 2 shows the locus where they provide the best match 
between the relic and local BHMFs. The solid and dashed lines shows the loci 
where the average deviation between the BHMFs is less than 1 and 0.7a, respec- 
tively. Outside of the solid contour, the agreement between the BHMFs is poor. 
The dotted lines marks the e values for a non-rotating Schwarzschild BH and a 
maximally rotating Kerr BH. Acceptable values are in the range e = 0.04 — 0.16 
and A = 0.1 - 1.7. 

The discrepancy found by [22,3] is thus removed without requiring large 
accretion efficiencies because of (i) the use of the zero points of [11], (ii) the 
non-zero intrinsic dispersion of the Mbh-c* and MeH-ibui relations, and (iii) 
the smaller bolometric corrections. While (ii) and (iii) are generally accepted, 
(i) is more controversial because the definition of 'secure' BH mass measurement 
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Fig. 3. (a) Average growth history of BHs. The symbols indicate the points when a BH 
reaches a given fraction of its final mass, (b) Average mean lifetimes of active BHs as a 
function of their mass at z = 0. The solid line corresponds the standard case (e = 0.1, 
A = l). 



varies from author to author. However, even using the same zero points as [22] 
to estimate the local BHMF we obtain e = 0.05 - 0.2 and A = 0.15 - 2.5, i.e. 
efficiency is still e < 0.2. It is worth noticing that the average limits on L/L-^dd 
that we find are perfectly in agreement with the average values estimated by [14] 
on a large sample of SDSS quasars (see their Fig. 2). 

5 Anti-hierarchical BH growth and AGN lifetimes 

Fig. 3a shows the average growth history of BHs with different starting masses 
at z — 3. Symbols mark the point when a BH reaches a given fraction of its final 
mass. At z < 3, all BHs gain at least 95% of their final mass but BHs which 
are more massive than 10 8 M Q grow earlier and gain 50% of their final mass 
by z ~ 2. Smaller BHs grow at lower redshifts (z < 1). This anti-hierarchical 
growth of BHs is a consequence of the redshift evolution of the [21] LF. 

Fig. 3b shows the average total lifetimes of active BHs, i.e. the time required 
for the BH growth since z = 3. The solid line shows the "canonical" case with 
e = 0.1, A = 1. Lines with symbols show limiting cases from Fig. 2. Local high 
mass BHs (Mbh < 10 9 M ) have been active, on average, ~ 1.5 x 10 8 yr. On 
the contrary, the assembly of lower mass BHs has required active phases lasting 
at least three times that much (~ 4.5 x 10 8 yr). The average lifetimes can be 
as large as 10 9 yr with the smaller e and A values compatible with local BHs 
(e = 0.04, A = 0.1 - see Fig. 2). 

Overall, the plots in Fig. 3 indicate that smaller BHs (Mbh < 10 s Mq) 
find more difficulties in growing than larger ones. Indeed, this is consistent with 
physical models for the coevolution of BHs and galaxies. Smaller BHs form in 
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shallower potential wells with respect to more massive ones and are thus more 
subject to feedback from star formation (e.g. supernovae explosions) and from 
the AGN itself [15]. 

6 Conclusions 

We have shown that the local BH mass function and that of AGN relics are in 
good agreement with standard accretion efficiency and L/L^dd ratio (e ~ 0.1, 
A ~ 1). In particular, the local BH Mass function implies that the density in 
BHs is pbh = 4.6(-1.4;+1.9)(/i/0.7) 2 x 10 5 M Q Mpc" 3 , a factor 1.8 higher 
than estimated by [22]. Merging of BHs is cither not important or it does not 
significantly alter the relic BHMF, at least at z < 3. The BH growth is anti- 
hierarchical, in the sense that smaller BHs (Mbh < 10 7 M©) grow at lower 
redshifts (z < 1) with respect to more massive ones (z = 1—3). The global picture 
which emerges is than that local BHs grew during AGN phases in which accreting 
matter was converted into radiation with e = 0.05 — 0.2 and L/Lboi = 0.15 — 2.5. 
The average total lifetime of these active phases ranges from ~ 4.5 x 10 8 yr 
(M BH < 10 7 M©) to ~ 1.5 x 10 s yr (M BH > 10 9 M ). 
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